P-wave Radial distributions of a Heavy-light meson on a lattice 
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Abstract 

This is a follow-up to our earlier work for the charge (vector) and matter (scalar) distributions 
for S-wave states in a heavy-light meson, where the heavy quark is static and the light quark has a 
mass about that of the strange quark. The calculation is again carried out with dynamical fermions 
on a 16 3 x 24 lattice with a lattice spacing of about 0.14 fm. It is shown that several features of 
the S- and P-wave distributions are in qualitative agreement with what one expects from a simple 
one-body Dirac equation interpretation. 
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I. INTRODUCTION 



Experimentalists are often able to tell us several properties of a given meson, such as its 
energy, width and angular momentum. However, usually they can not tell us the structure 
of the meson. For example, with B s states — the topic of interest here — they can not 
say whether these states are bs, bsuu or BK. Unfortunately, when theoreticians try to 
describe B s states, they have to decide beforehand the state structure to be used in some 
model, which often has sufficient freedom to fit the data with any of the possible structures. 

In an attempt to clarify the Experiment « ► Theory comparison, we suggest the use of 

lattice QCD. In principle, lattice QCD should give us all we need to know about B s states. 
However, in practice, the results need to be corrected for the lattice spacing (a), finite lattice 
size (L) and quark mass (m q ) effects — but these are usually under control and with the 
advent of more computer resources they will decrease in importance. 

The strategy followed here is to concentrate on the simplest of quark states, namely, the 
Qq system, where Q is an infinitely heavy quark (i.e. static) and q is a quark with about 
the strange quark mass. This system is sufficiently simple to enable state-of-the-art lattice 
calculations to generate much more "data" than can be achieved by direct experiment. This 
data consists of the ground state energies of S-, P- and D-wave states and also the spin 
average F-wave energy Not only are the ground state energies extracted, but also those 
of the corresponding excited states containing at least one radial node. In addition, the 
vector (charge) and scalar (matter) radial distributions of these states can be measured. 
This is an abundance of data, far beyond what has been done experimentally in the B s 
meson. Given all this data the challenge is now for theorists to make models to explain 
it. In this quest there are two simplifications. Firstly, since Q is static, the system is 
essentially reduced to a one-body problem involving only the light quark q s . Secondly, 
in the lattice calculations, since the energies and radial distributions are extracted from Qq 
correlations propagating in Euclidean time, it is expected that the resultant states are indeed 
Qq states with little contamination from other possible multiquark components. Support 
for this expectation (hope) is seen in Fig. 2 of Ref. Q|. There a QQ correlation generates 
a linearly rising potential for interquark distances far larger than expected i.e. way beyond 
where (Qq) (Qq) configurations should appear through string breaking. To see this effect the 
explicit introduction of (Qq)(Qq) correlations is needed 0. 



2 



The "data base" for properties of Qq states measured on lattices is so far incomplete. 
In Ref. (4 1 we concentrated on the S-wave energies and radial distributions. The latter 
showed two distinct features. Firstly, the excited state distributions exhibited nodes as 
expected. Secondly, the charge distribution (xc(R)) was of longer range than the matter 
distribution (xm{R)) but with xc(0) ~ x^(0). This is readily explained by the one-body 
Dirac equation, since there xc{R) = G(R) 2 + F(R) 2 and xm{R) = G(R) 2 — F(R) 2 , where for 
S-waves G(0) ^> P(0) with G(R) and F(R) becoming comparable for large R. In Ref. jlj the 
data base was extended to include the energies of the excited states P+/~ and D + /_ and also 
the spin averaged combinations P+_ and F+-. Here P is the P-wave state with j q = 1/2, 
since the spin of the Q does not play a dynamical role. In this note we return to measuring 
the ch arg e and .natter distributions as in Re, Q, but eoneentrate on the P_-state end its 

n 

excitations. This state is of particular interest, since — as shown in Fig. 4 of Ref. — the 
indications are that the predicted B s (0 + ) is below the Pi^-threshold and so should be very 
narrow as was found for the cs counterpart. The outcome is seen in Figs. 1 for the ground 
state distributions x 11 and the off-diagonal distributions between the ground state and the 
first excited state x 12 . 

The latter show single nodes as expected from a first excited state. However, x 11 shows 
two features which at first sight seem surprising: 

1) The distributions are finite at R = 0, even though they are P-waves. 

2) The matter distribution has a node, even though it involves only the ground state. 
But again this is precisely what one expects from solutions of the Dirac equation, where for 
the P_-state both G and F are non-zero at R = and, furthermore, the "small" component 
F can be larger than G at small R. In Figs. 2 the above data, now with a factor of R 2 
included, are compared with the Dirac distributions using a quark mass of 100 MeV and 
an interquark potential V = —a/R + bR, where a = 0.6 and 6=1.3 GeV/fm. These three 
parameters are very sensible and were not tuned to get an optimal fit. 



II. CONCLUSIONS 



The S- and P-wave charge and matter distributions in Figs. 2 suggest that they can be 
understood qualitatively in terms of the one-body Dirac equation. The challenge is now 
to see to what extent there is an analogous quantitative description of the energies and 
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FIG. 2: The S and P_ radial distributions from Fig. 1 compared with Dirac equation fits. 

distributions of all the S-, P-, D- and F-wave states — both ground and excited, when they 
become available. It is possible that the strategy used in studying the iVTV-potential is 
appropriate, namely, to first concentrate on the higher partial waves and so avoid or reduce 
complications, such as the effect of form factors needed to regulate the one-gluon-exchange 
potential and also instanton-induced interactions, that enter at small values of R. 
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